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The pulse tube engine is a simple heat engine based on the pulse tube process. Due to its simplicity it has 
a high potential to be applicable in waste heat usage and energy harvesting purposes. In this work, 
mathematical and experimental design tools are developed to study a pressurized laboratory scale pulse 
tube engine. The mathematical model is based on the transient numerical solution of the governing 
differential equations for mass, momentum and energy. The Modelica environment of SimulationX is 
used to solve the equations numerically and the model is employed to design the experimental test 
engine with helium as working fluid. The transient behavior of the pulse tube engine's underlying 
thermodynamic properties is studied numerically and experimentally under different design parameters 
as well as for different heat input temperatures, filling pressures and operating frequencies. The 
measured engine characteristics are compared with the calculated predictions. Internal and external 
power losses are quantified. Design studies for a further development of the pulse tube engine are 
performed experimentally. The developed numerical tool provides a rational framework for up-scaling 
the current laboratory model to industrial scale. 

© 2013 Elsevier Ltd. All rights reserved. 


1. Introduction 

The pulse tube engine is a simple prime mover for the conver¬ 
sion of externally supplied thermal power into mechanical power 
and subsequently into electrical power by using a generator. It 
represents the thermodynamic counterpart of the pulse tube 
refrigerator process presented in the 1960’s by Gifford and Longs- 
worth [1], Similar engines were described by Chen at al. [2] and 
Tailer [3], The concept of the pulse tube engine was presented by 
Organ [4] and Hamaguchi et al. [5]. Several experimental in¬ 
vestigations to characterize this engine were also published by 
Hamaguchi et al. [5—7]. A theoretical analysis of pulse tube engine 
was done by Moldenhauer et al. [8]. It was claimed that a broken 
thermodynamic symmetry is the underlying working principle of 
the pulse tube engine. Thermodynamic asymmetry results from the 
fluidic linkage of an adiabatic space (pulse tube) and an isothermal 
space (regenerator). 

The functional principle of the pulse tube engine is based on the 
simultaneous compression and motion of the working gas into an 
area of high temperature of the wall material, followed by a 
simultaneous expansion and motion of the working gas into an area 
of low temperature of the wall material. 
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In spite of the recent modeling attempts [8], the design of pulse 
tube engines is still largely a matter of intuition and its working 
mechanisms are not sufficiently well understood. In particular, it 
would be desirable to possess an experimentally validated simu¬ 
lation tool that would allow reliable extrapolation of the perfor¬ 
mance of pulse tube engines to higher power ratings. The goal of 
the present paper is to provide experimental and numerical tools to 
characterize the mode of operation of the pulse tube engine under 
different heat input temperatures, filling pressures and operating 
frequencies, and to find optimal design parameters for future high- 
performance pulse tube engines. 

In Section 2 the functional principle of the pulse tube engine is 
presented. The developed numerical simulation model is given in 
Section 3 and the design of the experimental pulse tube engine in 
Section 4. After the model calibration in Section 5, the experimental 
test environment is characterized in Section 6. A comparison of the 
measured and calculated engine performance is presented in Sec¬ 
tion 7 and the observed influence of design parameters is given in 
Section 8. 

2. Basic design and functional principle of the pulse tube 
engine 

The basic design of the investigated pulse tube engine is shown 
in Fig. 1. A compression cylinder supplies a pressure swing and gas 
shifting to the connected system consisting of an adiabatic buffer 
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List of symbols 

Lr 

piston rod length [m] 



Lpt 

pulse tube length [m] 

a 

heat transfer coeff. [W/(m 2 K)] 

Lreg 

regenerator length [m] 

At 

time increment [s] 

m 

gas mass [kg] 

E 

energy flux [W] 

M br 

brake momentum [N m] 

Ecyc 

cycle energy flux [W] 

Mdr 

drive momentum [N m] 

H 

enthalpy flux [W] 

n 

number of subvolumes 

m 

mass flux [kg/s] 

n p 

polytropie exponent 

Q 

heat flux [W] 

Nu 

Nusselt number 

Qcg 

cond. heat flux gas [W] 

P 

pressure [Pa] 

Qcm 

cond. heat flux mat. [W] 

Pf 

filling pressure [Pa] 

Q-cyc 

cycle heat flux [W] 

Pg 

gross power [W] 

Qext 

external heat flux [W] 

Pn 

net power [W] 

Q-hx 

gas-material heat flux [W] 

Pi,ext 

external power loss [W] 

Q-in 

supplied thermal power [W] 

Pl.int 

internal power loss [W] 

W 

pV — power [W] 

Pmax 

maximum pressure [Pa] 

Vn 

net efficiency 

Pmin 

minimum pressure [Pa] 

Vsh 

shaft efficiency 

Psh 

shaft power [W] 

K 

adiabatic exponent 

Pr 

Prandtl number 

X 

heat cond. gas [W/(mI<)[ 

r c 

compression ratio 

p 

dynamical viscosity [Pa s] 

r d 

distribution ratio 

GJ 

cycle frequency [Hz] 

R s 

spec, gas constant [J/(kg K)] 


piston starting position [deg] 

Ecr 

crank ratio 

n 

pressure ratio 

Re 

Reynolds number 


heat transfer module 

s 

stroke length [m] 

p 

gas density [kg/m 3 ] 

T 

gas temperature [°C] 

Pm 

material density [kg/m 3 ] 

t 

time variable [s] 

T 

temperature ratio 

T c 

cold temperature [°C] 

TO 

minimal temperature ratio 

T h 

hot temperature [°C] 

A, B 

fitting parameters [Ws, Ws 3 ] 

Tm 

material temperature [°C] 

a, b 

fitting functions [Wj 

Tchx 

gas temperature chx [°C] 

Af 

free flow area [m 2 ] 

Tcyl 

gas temperature cylinder [°C] 

Am 

material area [m 2 ] 

Thhx 

gas temperature hhx [°C] 

Ahx 

heat transfer area [m 2 ] 

Tpr 

gas temp, pulse tube [°C] 

Cm 

spec, heat capacity [J/(kg 1<)[ 

Treg 

gas temp, regenerator [°C] 

C P 

spec, heat capacity [J/(kg 1<)[ 

u 

velocity [m/s] 

c v 

spec, heat capacity [J/(kg 1<)[ 

u m 

internal energy material [J] 

d 

regenerator diameter [m] 

V 

gas volume [m 3 ] 

d h 

hydraulic diameter [m] 

V 0 

dead volume [m 3 ] 

d w 

wire diameter [m] 

v d 

displacement volume [m 3 ] 

e 

porosity 

V m 

material volume [m 3 ] 

F 

force [N] 

V hx 

gas vol. heat exchangers [m 3 ] 

f 

operating frequency [Hz] 

*max 

max. engine gas volume [m 3 ] 

ff 

flow friction factor 

Vmin 

min. engine gas volume [m 3 ] 

F f 

flow friction force [N] 

v pc 

gas volume pulse tube [m 3 ] 

I 

momentum [N m] 

Vreg 

gas volume regenerator [m 3 ] 

L 

control volume length [m] 

X 

space variable [m] 


space (pulse tube) and an isothermalizer (regenerator). Both com¬ 
ponents are fluidically linked via a gas heater (hhx: hot heat 
exchanger) at the temperature T b , providing the heat input. A gas 
cooler fcliX]: cold heat exchanger 1) at the temperature T c is located 



Fig. 1. Basic design of the pulse tube engine: chxi - cold heat exchanger 1, chx 2 - cold 
heat exchanger 2, hhx - hot heat exchanger. 


between cylinder and pulse tube. To set the depicted material 
temperature distribution, a heat sink (chx 2 i cold heat exchanger 2) 
at the temperature T c is located at the end of the regenerator. 

The conversion of thermal power into mechanical power is 
enabled by the simultaneous compression and motion of the 
working gas into an area of high material temperature, followed by 
a simultaneous expansion and motion of the working gas into an 
area of low material temperature. Thus, the mean pressure during 
expansion is higher than during compression resulting in a release 
of mechanical power. In Ref. [8] was shown that the pulse tube 
engine is working as a free piston engine. However, enabling the 
usage of a simple electrical generator, the linear motion of the 
piston is converted into rotation of a flywheel. 
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3. Numerical simulation model 

A control volume based simulation model is developed. The 
model is an extension of the ideal pulse tube engine model pre¬ 
sented in Ref. [8], It takes into account dissipative mechanisms such 
as friction, heat shortcuts due to heat conduction and non-ideal 
heat transfer. 

Fig. 2 shows the schematic representation of the pulse tube 
engine. Heat can be exchanged between the surroundings and the 
engine material via the heat flux Q ext as well as between the engine 
material and the working gas via the heat flux Q hx . In the cylinder 
the working gas is charged with the compression/expansion power 
W. The components are connected at equal pressure through mass 
fluxes m and enthalpy fluxes H. Due to the high temperature gra¬ 
dients, the pulse tube and the regenerator are each divided into sub 
control volumes (one-dimensional approximation). The other 
components consist of one control volume each (zero-dimensional 
approximation). All entering fluxes are given a positive sign. The 
following simplifications are made: 

1) The flow is one-dimensional. 

2) The working fluid is a perfect gas with p = pR s T, constant spe¬ 
cific heat capacities c p , c v and the specific gas constant R s . 

3) The thermodynamic and fluid mechanical properties are uni¬ 
form within a given control volume. 

4) The speed of sound in the working fluid helium has in the 
temperature range between 20 °C and 600 °C a value higher 
than 1000 m/s. The typical operating frequency of the pulse 
tube engine is lower than 100 Hz. For a piston stroke of 50 mm 
the maximum piston velocity is thus about 10 m/s. The velocity 
of the fluid has the same order of magnitude. Since the speed of 
sound is two orders of magnitude higher than the velocity of 
the fluid, pressure wave reflections at system boundaries are 
neglected. 


With the constant material density p m and the constant specific 
heat capacity of the material c m the differential equation for the 
material temperature T m can be written in the form 

dT' 

PmVmCm-^jr = “Qhx + Qcm + Qext- (4) 

Eqs. (2) and (4) contain additional diffusive terms, Q cg and Q cm , 
describing the heat conduction within the working gas and the 
material, respectively. 

The heat flux Q hx between working gas and material is calcu¬ 
lated with the heat transfer area Ah x through 

Qhx = a ^hx(Jm — T), (5) 


where the heat transfer coefficient is calculated using a = NuX/dh, 
with the Nusselt number Nu, the heat conductivity of the gas X and 
the hydraulic diameter d/,. The Nusselt number is obtained from 
experimental correlations depending on the flow field. For the 
cylinder, the heat exchangers and the pulse tube, the heat transfer 
coefficient a is obtained from experimental correlations given in 
Ref. [9], Since the porous structure of the regenerator is made of 
stacked wire mesh screens, experimental correlations for a from 
Ref. [10] are used. The heat transfer area Ah X in the regenerator with 
the length L reg and the diameter d reg is given by 


A hx = e(l -e)- 


nd?~ a L, 


reg L reg 


( 6 ) 


where e is the regenerator porosity and d w the wire diameter of the 
mesh. The heat transfer coefficient is calculated according to [10] 
from 


a 


c p p 
d h Pr 2 / 3 ’ 


(7) 


3.1. Basic equations and gross power output 


The differential equations for pressure p and temperature T of 
each control volume V are obtained from the ideal model in Ref. [8] 
with 


dp 

dt 


RJ 

V 


. p dV pdT 
I>-l7dF + 7=dF 


( 1 ) 


= -^Em + Qhx + Qcg + EH + W. (2) 

For the material of the control volume V m , the conservation 
equation of the material's internal energy U m is given by 

= ]TQ = “Qhx + Qcm + Qext- (3) 



with the dynamical viscosity p, the Reynolds number Re and the 
Prandtl number Pr. The heat transfer module \j/ is calculated ac¬ 
cording to 

i = (2.128- 8.591e + 9.6878e 2 )Re-°' 113 -° 475e . (8) 

The cylinder volume varies in time and is given by 

V(t) = V 0 + y ^1 + r cr - cos(wt - cp) - \J r 2 r - sin 2 (wt - <p)\, 

O) 

where the crank ratio r cr = 2 L r /s is the ratio between piston rod 
length L r and crank radius s/2, which is half of the stroke s. The cycle 
frequency is defined by to = 2te f and the starting position of the 
piston is given by <p, where <p = 0 at the piston’s top dead center 
(TDC). The piston’s bottom dead center is denoted with BDC. The 
dead volume and the displacement volume are represented by Vo 
and Vd. 


CV 
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Fig. 2. Schematic representation of the pulse tube engine divided into control volumes 
with mass and energy fluxes crossing the control volume borders. 


Fig. 3. Acting pressure forces, momentum forces and the friction force on the control 
volume gas mass. 1 = left side, r = right side. 
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Using Eqs. (1)—(4), the gross power output P g of the pulse tube 
engine can be calculated. 

3.2. Net and shaft power output 

The engine’s gross power output is reduced by internal power 
losses P/.int, mainly due to friction between gas and filling ma¬ 
terial as well as heat transfer. In order to calculate the net power 
output P n of the engine, the equation of motion is added to the 
model. 

The momentum balance is given by 

l-Ef. no) 

where I = mu represents the momentum of the gas element m 
having the velocity u. Forces on the gas element are represented by 
F. With the expression m = pAfU, which relates the mass flux m 
with the velocity, the free flow area Af and the gas density p = m/V 
follows 


To get the shaft power output P s h, the net power has to be 
reduced by an amount equaling the power losses due to friction in 
the cylinder, the bearings and the crank mechanism. The friction 
force in the cylinder is proportional to the constant gravity force of 
the piston and piston rod. For the crank mechanism the friction 
force in the linear and rotating bearings is a superposition of the 
constant gravity force of the rod and the shear force caused by the 
crank mechanism. The shear force is proportional to the square of 
the operating frequency/ 2 . Thus, the total friction power loss P/ >ex t 
in the cylinder and the crank shaft can be modeled according to 

P/,ext(f) = Af + Bf 3 , (16) 

where the influence of the filling pressure variation is neglected 
and A and B are constant fitting parameters. The shaft power output 
can then be calculated from P s h = P n + Pi,ex t. with the convention 
that any power delivered by the engine has a negative sign. 

The model has two free parameters which need to be quantified. 
These are the dead volume V 0 and the starting position of the piston 
<p in degree of crank angle relative to TDC. Both parameters will be 
determined experimentally. 


m Vm rh 

pAf mAf m' 


( 11 ) 


where L represents length of the control volume (cv). Using Eq. (11 ), 
the momentum results in / = mu = rhL and thus Eq. (10) in 


d/ .dm ^ 

di =L W = 2^ F 


( 12 ) 


According to Fig. 3 momentum forces are calculated with mu = 
Lift 2 /m and pressure forces with pAf. The flow friction force is 
represented by F/. Hence, together with the thermal state equation 
of a perfect gas, Eq. (12) becomes 


dm 

dt 


■ iRsTt ■ 2 RsTr Af 

m ‘W~ mr W + T {Pl ~ Pr) 



(13) 


According to [10] the flow friction force can be modeled using 


F f 


= ff 


2rh|rh|L 

pAfd h 


(14) 


3.3. Computation parameters 

The system of differential equations is solved numerically using 
the object-oriented programming language Modelica within the 
software SimulationX. With the correct initial conditions, the 
dynamical state variables of the pulse tube engine such as the gas 
temperatures {T cy i(f), T chX] (t), T chx? (t), T hhx (t), T pt (x,f), T reg (x,t)} are 
obtained. The cycle energy fluxes £ C yc are calculated by integrating 
the energy fluxes E over one engine cycle and multiplying the result 
with the operating frequency/. As an example, the cycle heat flux of 
a component and the net cycle power are given by Q cyc = ff Qdt 
and P n = ff Wdt. For all calculations a sufficient number of engine 
cycles is calculated until the changes from cycle to cycle are 
negligible. The values of Q cyc and P n are averaged over the last ten 
cycles. 

The calculations are performed with minimum and maximum 
permissible time increments of At m j n = 1CT 8 s and Af max = 1(T 5 s, 
and a number of subvolumes for pulse tube and regenerator of 
n = ipt/reg/2, where L pt / r eg ate the lengths of the components 
measured in millimeters. 


where ff represents a flow friction factor. 

The temporal change of the mass flux m becomes finally 


dm 

dt 


R s fmfT, 

v{p, 




(Pi ~ Pr) 




2rh|rh| 

pAfd h 


(15) 


The use of the absolute value for the mass flux ensures that the 
direction of the flow friction force is opposite to the direction of 
motion. 

The gas flow in the cylinder, the heat exchangers and the pulse 
tube has low Reynolds numbers (Re < 1000), so that the friction 
factor relation ff = 16/Re for laminar flow is applicable. In the 
regenerator even lower Reynolds numbers (Re < 100) are present. 
According to [11,12] the flow field in porous media can be turbulent 
at low Reynolds numbers, too. To calculate the friction force in the 
regenerator an experimental friction factor correlation from Ref. 
[10] is used. The friction factor relation ff = 16/Re and the friction 
factor correlation from Ref. [10] are valid for fully developed sta¬ 
tionary flows, but also used for the present transient analysis. The 
induced modeling error resulting from this is judged to have minor 
importance compared to the error induced by the one-dimensional 
modeling. 


4. Design of the experimental pulse tube engine 

The simulation model is used to design an experimental test 
engine. The aim of the experimental set-up is to validate the 
simulation model. Furthermore, the experimental test engine is 
also used to study the pulse tube engine’s performance for different 
designs and operation conditions. 

4.1. General engine design 

The general design of the test engine is shown in Fig. 4. It has a 
cylinder with a bore of 24 mm and a variable stroke between 
20 mm and 100 mm. The cylinder is a low friction A1RPEL (Airpot 
cooperation) pneumatic cylinder with a glass liner and a piston 
made of aluminum and graphite. The front side of the cylinder is 
connected to a water cooled copper heat exchanger (chx!). The 
pulse tube, the hot heat exchanger (hhx) and the regenerator are 
placed in a stainless steel tube with an inner diameter of 20 mm 
and a total length of 230 mm (including flanges). An electrical 
resistance heater is place inside the hot heat exchanger. At the cold 
side of the regenerator a second water cooled heat exchanger (chx 2 ) 
is located. 
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4.2. Model based design of the components 

v max = V d + V 0 + V p t + Vreg + V hx 

Two parameters characterizing the design of the pulse tube 
engine — the compression ratio r c and the distribution ratio r d — are 

introduced. The compression ratio is defined as y min = Vq + V pt + \7 reg + \/ hx , 


(19) 


( 20 ) 


. _ ^pt + Weg + V'hx + Vd _ -i _ _ Yd _ 

Vpt + V'reg + V hx Vpt + V re g + V hx 

and the distribution ratio as 


(17) 


U- = 


Vreg 

Vpt ’ 


(18) 


where V r eg and V pt are the regenerator and the pulse tube volume, 
V d is the displacement volume of the cylinder and Vi lx is the sum of 
the heat exchanger volumes. Engine power output and efficiency 
have strong dependencies on the optimal choice of r c and r d . The 
present engine was designed to vary these values in a later opti¬ 
mization process. The distribution ratio was first set to r d = 1.51. 
Using a stroke of 40 mm the compression ratio is r c = 0.38. 

The regenerator is located at one side of the stainless steel tube. 
Therefore it is filled with stacked stainless steel wire mesh screens 
having a mesh number of 50 and a wire diameter of 200 pm. The 
length of the regenerator is 122 mm. Its porosity is 0.67 and the 
total regenerator mass is about 200 g. In the middle of the tube the 
hot heat exchanger with a length of 40 mm is located. The rest of 
the tube (68 mm) is used as pulse tube. A spacer with 10 mm in 
diameter and 40 mm in length is placed on the pulse tube axis in 
order to reduce its gas volume. Both water cooled cold heat ex¬ 
changers are made from copper. The cold heat exchanger chx 2 at 
the regenerator’s cold side has a 2 mm in diameter bore and a 
closed end. The cold heat exchanger chxi with a length of 11 mm at 
the pulse tube’s cold side has a spark-eroded structure to achieve a 
certain heat transfer structure. 

Heat input is provided by an electrical resistance heater and 
delivered to the working gas via the hot heat exchanger. Its heat 
transfer area has no theoretical limitation and was designed to be 
as large as possible. In contrast, the cold heat exchanger chxi needs 
a certain heat transfer area to allow optimal engine operation. To 
get the optimal heat transfer area, which maximizes the engine’s 
power output, the simulation model was used. The calculated 
optimal heat transfer area is 56 cm 2 . 

The resulting values of the gas and material volume, V and V m , as 
well as of the free flow area Ap the material area A m , the heat 
transfer area Ah X and the hydraulic diameter d d are listed in Table 1. 
The maximum and minimum gas volume of the engine are given by 


where V 0 represents the sum of the dead volumes of the cylinder as 
well as of the connected tubes and valves. The dead volume Vo is at 
low temperature and can entirely be regarded as cylinder dead 
volume. It is split into a known (measurable) dead volume Vg and 
an unknown dead volume Vg. The measured known dead volume is 
Vg = 7.5 cm 3 and the free parameter V 0 is replaced by the free 
parameter Vg. Its value is calculated from pressure measurements 
after building up the engine. According to [ 8 ] it partly represents 
virtual dead volume due to leakage between piston and cylinder 
liner. 


4.3. Experimental set-up and operating conditions 

The structure of the experimental set-up is shown in Fig. 4. The 
rear side of the cylinder is connected to an external gas volume 
simulating a pressurized environment. The engine is filled with 
pressurized helium up to a filling pressure ofp/= 12 bar. The piston 
front and rear side have a connecting tube, which can be opened 
and closed using valve Vp This bypass tube allows to run the cyl¬ 
inder without any pressure swing in the system and almost no gas 
passing chxp The cold heat exchanger chx 2 can be opened and 
closed to the gas volume using valve V 2 . This allows to fill the en¬ 
gine with pressurized gas without having Vi opened. The whole 
system is connected to a vacuum pump via valve V 3 and a com¬ 
pressed gas cylinder via valve V 4 . The system can be emptied to the 
ambient by opening valve V 5 . The filling pressure is measured using 
a mechanical pressure gauge with a precision of 0.25 bar. 

The piston rod is guided by linear bearings. The linear move¬ 
ment is transformed into rotation via a crank shaft mechanism. The 


Table 1 

Geometrical parameters of the engine components. 


Parameter 

chx 2 

Regenerator 

hhx 

Pulse tube 

chx! 

Cylinder 

V [cm 3 ] 

1.8 

25.6 

2.8 

16.9 

0.8 

21...39“ 

Vrr, [cm 3 ] 

- 

20.8 

2.4 

4.5 

1.8 

- 

Af [cm 2 ] 

0.1 

2.1 

0.7 

2.4 

0.8 

4.5 

A m [cm 2 ] 

- 

1.7 

2.4 

0.7 

1.6 

- 

Ahx [cm 2 ] 

18 

1700 

187 

43 

56 

44...74 a 

dh [mm] 

4 

0.4 

0.6 

10 

0.6 

24 


a Oscillating. 
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connected flywheel serves as rotational energy storage for the 
compression stroke. Two electrical motors are connected to the 
flywheel. One is used as electrical dynamometer to apply a brake 
momentum on the engine. The resulting electrical current is 
measured and used to calculate the brake momentum Mt, r . The 
second electrical motor serves as a drive to motor the engine. The 
electrical current of the drive motor is also measured to calculate 
the drive momentum Ma r . 

The electrical heater, placed inside the hot heat exchanger has a 
maximum permissible temperature of 700 °C and a maximum 
thermal power of 160 W. The heater temperature is measured and 
controlled using a two-level controller. It can be fixed between 
40 °C and 650 °C. The cold heat exchangers are water cooled with a 
cooling water inlet temperature of 20 °C. 

4.4. Measurement techniques 

In order to characterize the pulse tube engine in operation, the 
signals of various sensors and measuring devices are monitored. 
The main measured quantities are temperatures T, the supplied 
thermal power Q in , brake and drive torques JVfb r and JV?d r , the 
working gas pressure p, the operating frequency / and the crank 
angle of the flywheel ip. The measured temperatures are in partic¬ 
ular the temperatures of the heat exchangers T chX] , T chX2 and Thhx. as 
well as the temperature profiles of the regenerator Treg and the 
pulse tube Tpt.A scheme of the measurement is shown in Fig. 5. 

The thermal power input Q in is calculated from the measured 
voltage and current of the heater power supply. Heat losses are 
quantified experimentally in Section 6.1. The torques of the two 
electrical motors Mb r and Md r are calculated from its electrical 
currents, which are correlated to measured voltages of shunts. The 
current of the brake motor is smaller than of the drive motor. To 
increase the accuracy three switchable shunts are used. The 
measured voltages are amplified and processed by a micro 
controller. To ensure the highest possible precision, the amplifica¬ 
tion circuits were checked for their linearity before the measure¬ 
ments. The conversion factors between the measured voltages and 


the quantities to be identified were determined by calibrations. The 
values are stored in the micro controller Flash-EEPROM. 

The working gas pressure is recorded using piezoresistive 
pressure sensors p i and p 2 between chxi and the cylinder as well as 
at chx 2 . The sampling rate is 4 kHz. The voltage signals are logged 
by a data acquisition system and monitored with a personal com¬ 
puter (PC). Operating frequency /and crank angle <p are measured 
by two light barriers at the interrupt input of the micro controller. 
The first light barrier is interrupted 20 times per revolution of the 
motor shaft, thus permitting a frequency measurement with an 
accuracy of 0.05 Hz. The determination of the crank angle becomes 
possible with the second light barrier which is interrupted only 
once per revolution and thereby setting the zero position of the 
crank angle. Each other interruption of the light beam of the first 
barrier increases the value of the crank angle by 18°. The crank 
angle is converted into a voltage between zero and 5 V via pulse 
width modulation and then forwarded to the data acquisition 
system and the personal computer. With this data the pV-diagram 
of the engine can be monitored. 

The temperatures are measured with K-type thermocouples. 
The thermoelectric voltages are converted into temperatures by an 
AGILENT 34970A transducer with built-in reference junction. The 
data are transmitted via the RS232 interface every 2 s to the per¬ 
sonal computer. The thermoelectric voltage of the thermocouple 
placed at the hot heat exchanger is also amplified externally and 
read with the micro controller. The signal is used to control the 
heater temperature with a two-level controller. 

The errors of the measured quantities are listed in Table 2. For 
the non-derived quantities/, pp, T c h x , Thhx. T pt , T reg , M b r , M<j r and Q in 
the listed relative or absolute errors are systematic errors of the 
measurement set-up. For the derived quantities P S h, P n , Vsh and p n 
the systematic errors are obtained through the law of error prop¬ 
agation. Using the electrical dynamometer the minimum absolute 
error of the shaft power is quantified experimentally to 0.1 W. 
Measuring the shaft power or the net power with the electrical 
drive motor the minimum absolute error is quantified experi¬ 
mentally to 0.4 W. 



Fig. 5. Structure of the data acquisition system. 
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Table 2 

Systematic errors of the measured quantities. 


Measured quantity 

Max. relative error 

Min./fix absolute error 

/ 

1% 

Min. 0.05 Hz 

Pf 

- 

Fix 0.25 bar 

Tchx 

- 

Fix 1 K 

Thhx 

- 

Fix 15 K 

Tpt, Treg 

- 

Fix 30 K 

Mbr. Mdr 

10% 

— 

Qin 

4% 

- 

Psh 

11% 

Min. 0.1 W a 

Psh, Pn 

11% 

Min. 0.4 W b 

Vsht Vn 

15% 

- 


a Using the electrical dynamometer. 
b Using the electrical drive motor. 


5. Model calibration 

The free parameters which need to be quantified are the un¬ 
known dead volume Vfi (see Section 4.2) and the starting position 
of the piston q> (see Section 3). Both values were obtained from the 
measured time-dependent pressure at ambient temperature. The 
parameters Vft and rp were varied to minimize the error of the 
calculated minimum pressure and the calculated pressure ampli¬ 
tude for filling pressures between 6 bar and 10 bar and a mean 
engine frequency of 8 Hz. The results are <p = 76 deg and 
V$ = 13 .5 cm 3 . Thus, the total additional dead volume of the en¬ 
gine is Vo = 21 cm 3 . A comparison of the calculated and measured 
time-dependent pressure is shown in Fig. 6. The relative differences 
between the measured and calculated minimum pressure and 
pressure amplitude are smaller than 2.5% and 3.2%, respectively. 

6. Characterization of the experimental test environment 

To characterize the experimental set-up and to quantify the 
measured power values, the thermal power loss of the electrical 
heater and the power dissipation in the cylinder and the crank 
mechanism are investigated. 

6.3. Measured versus calculated thermal power input 

The engine was run at 7), = 600 °C, T c = 25 °C and filling pres¬ 
sures between 6 bar and 10 bar. The electrical dynamometer was 
used to apply a brake momentum on the engine. A comparison of 
the measured and calculated thermal power input is given in Fig. 7. 
The power input was measured from the electrical power con¬ 
sumption of the heater placed inside the hot heat exchanger. The 
measured thermal power input is between 40% (10 bar, 8 Hz) and 
95% (6 bar, 4 Hz) higher than the calculated thermal power input. 
The reason is that not all the heat power is transferred to the 
working gas. A part of it is lost to the ambient air. However, 
calculated and measured heat power input increase with the 
operating frequency at an average slope of about 7.2 W/Hz. 

The calculation predicts a monotonic increase of the thermal 
power input with the filling pressure. In the experiment the ther¬ 
mal power input increases only between 6 bar and 8 bar of filling 
pressure. For filling pressures higher than 8 bar, the measured 
thermal power input decreases. A reason for this behavior could be 
a smaller temperature difference between the working gas and the 
heat exchanger material and thus a lower entropy production 
during the heat transfer process. 

6.2. Friction power loss in the cylinder and the crank mechanism 

The friction power loss in the cylinder and the crank mechanism 
was measured for filling pressures between 6 bar and 10 bar and 



Fig. 6. Calculated and measured time-dependent pressure as functions of the crank 
angle for different filling pressures. 

frequencies between 4 Hz and 12 Hz. To do so, the engine was 
motored using the drive motor while valve V] was open. Thus, no 
compression occurred in the cylinder. The result is shown in Fig. 8. 
The observed influence of the filling pressure on the friction power 
loss is small. According to Section 3.2 the friction power was 
averaged over all filling pressures at each frequency and fitted using 
Eq. (16). The resulting fitting parameters are A = 0.301 ± 0.022 Ws 
and B = 0.005 ± 0.001 Ws 3 . 

6.3. Sign convention 

For all considerations in Sections 7 and 8 heat power inputs, 
mechanical power outputs and power losses have a positive sign. 
Mechanical power input by the drive motor is given a negative sign. 

7. Comparison of the measured and calculated engine 
performance 

The performance of the developed pulse tube engine was 
measured for different filling pressures, heater temperatures and 
operating frequencies. The observations are compared with the 
corresponding calculated predictions. To allow the identification of 



Frequency [Hz] 

Fig. 7. Calculated and measured thermal power input as functions of the operating 
frequency for different filling pressures. 
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Frequency [Hz] 

Fig. 8. Measured power consumption due to mechanical friction in the cylinder, 
bearings and crank mechanism as functions of the operating frequency for different 
filling pressures. To avoid too many lines in the plot the error bars are only shown for 
the measured power consumption at 8 bar. 

loss mechanisms all measurements were done for opened and 
closed valve V\. 

7.1. Measured shaft power output and efficiency 

The operating conditions of the engine were the same as in 
Section 6.1. The measured shaft power P S h and efficiency p s h for 
Th = 600 °C and filling pressures of 6 bar, 8 bar and 10 bar are shown 
in Fig. 9. The efficiency is calculated with i; sh = P S h/Qin> where Q in 
is the measured electrical power consumption of the heater. As 
shown in Section 6.1 the calculated thermal power input is between 
40% and 95% lower than the measured value. Thus, the thermal 
power delivered to the working gas is smaller leading to a lower 
efficiency than what was predicted. However, power and efficiency 
increase with increasing filling pressure. For a given pressure the 
peak frequency of the power is slightly higher than the peak fre¬ 
quency of the efficiency. Both values increase with increasing filling 
pressure. Having a filling pressure of 10 bar, the highest shaft power 
of about 0.9 W and an efficiency of about 1.2% is produced at a 
frequency of 5 Hz. According to Section 6.2 the friction power loss at 
a frequency of 5 Hz is 2.2 W. 



Fig. 9. Measured shaft power output and efficiency as functions of the operating 
frequency for different filling pressures. To avoid too many lines in the plot the error 
bars of the measured shaft efficiencies are not shown. 


7.2. Measured versus calculated net power output and efficiency 

The net power of the engine was obtained experimentally from 
the difference of the measured power consumption of the motored 
uncompressed engine (Vi opened) and the motored compressed 
engine (V] closed). In the uncompressed case, the drive motor has 
to deliver power P S h to the shaft, which equals the power loss due to 
mechanical friction P/ e xt in the crank shaft and the cylinder. In the 
compressed case, the amount of power delivered by the drive 
motor is reduced by the net power output P n produced by the pulse 
tube engine. The net power is calculated from P n = P s h + P/.ext and 
the net efficiency is calculated with i] n = P n /Q in , where P n is the 
measured net power and Q in the heat power input to the gas ob¬ 
tained from the calculation. Measurements and calculations were 
performed for heat input temperatures of 200 °C, 400 °C and 
600 °C, and filling pressures of 6 bar, 8 bar and 10 bar. 

A comparison of the measured and calculated results is given in 
Fig. 10. The filled symbols represent the measured and the unfilled 
symbols the calculated net power outputs and efficiencies. Within 
the restricted accuracy of a one-dimensional model, the calculated 
power outputs match well with the experimental results. At all 
filling pressures the relative difference in the calculated and 
measured net power outputs is for Th = 400 °C and Th = 600 °C 
smaller than 40%. Having a heat input temperature of Th = 200 °C, 
the relative difference is very high. A maximum error of 465% oc¬ 
curs at Th = 200 °C, Pf= 6 bar and/= 8 Hz. Due to its low absolute 
value, the uncertainty of the calculated net efficiency is quite high. 

The modeling error was estimated to about 20%. The maximum 
relative uncertainty of the measured net power is about 11 %, and 
since the electrical drive motor is used for the measurements the 
minimum absolute error of the measured net power is 0.4 W 
(compare to Table 2). The total relative difference of the net power 
outputs of calculation and experiment is thus about 31%. The error 
for Th = 400 °C and Th = 600 °C has the magnitude of this uncer¬ 
tainty. Since the net power output for T/, = 200 °C is smaller than 
0.4 W, the obtained relative error can be very high, but is still within 
the range of the measurement uncertainty. With increasing fre¬ 
quencies, the measured net power increases more strongly than the 
calculated net power. This could be due to an increase in the heat 
transfer rate because of flow transition from laminar to turbulent, 
which is not captured sufficiently by the model. 

The engine provides a net power output for all heat input 
temperatures and filling pressures, which increases with temper¬ 
ature and pressure. The increase slows down strongly for higher 
temperatures, such that the temperature has no significant influ¬ 
ence on the efficiency anymore. This result implies that with 
increasing temperature difference, more heat energy is dumped 
due to an enthalpy flow from the hot to the cold side of the pulse 
tube. For Th = 600 °C, p/ = 6 bar and / = 9 Hz the engine has a 
maximum net efficiency of about 8 %. The peak of the net power is 
always at a higher frequency than the peak of the net efficiency. 
Compared to the shaft power and efficiency, the net power and 
efficiency have their peak values about 5 Hz higher. This is due to 
the power dissipation in the cylinder and the crank mechanism (see 
Section 6.2). 

7.3. Internal power loss and gross power output 

The functional principle of the pulse tube engine is based on an 
intrinsically irreversible thermodynamic cycle [13]. Gross power is 
produced due to an irreversible thermal relaxation in the pulse tube. 

Fig. 11 shows the calculated pV-diagram of the pulse tube engine 
for pf = 10 bar, / = 8 Hz, Th = 600 °C and T c = 25 °C. The engine 
operates between minimum and maximum pressure, p m j n and 
pmax. respectively. Applying a global polytropic model [14] to the 
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(a) Net power output at Pf = 6 bar (b) Net efficiency at Pf = 6 bar 




(c) Net power output at pj = 8 bar 


(d) Net efficiency at Pf = 8 bar 




(e) Net power output at pf = lObar 


(f) Net efficiency at Pf — lObar 


Fig. 10. Calculated and measured net power output and efficiency as functions of the operating frequency for different filling pressures and heat input temperatures. To avoid too 
many lines in the plots the error bars of the measured net efficiencies are not shown. 


entire engine, a polytropic compression with the polytropic expo¬ 
nent iip(V) takes place from BDC toTDC. In the opposite direction a 
polytropic expansion with the polytropic exponent n e p (V) occurs. 
Since the function n^(V) differs from the function n£(V) the engine 
is able to provide a gross power output. 

With the ratio of the maximum to the minimum engine gas 
volume V'max/Vmin = 1-26 and the observed pressure ratio of 


If = 1.34, the mean polytropic exponent n p can be obtained 
from 

n= (v^y p (2i) 

V^min/ 

to be tip = 1.27. The working gas helium has a mean adiabatic 
exponent of k = 1.67. Thus, the global compression and expansion 
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processes tend to be close to isothermal behavior ( n p =1). Using 
this result, the gross power can be estimated with an isothermal 
model of the compression and expansion processes in combination 
with thermal (isochoric) relaxation occurring at TDC and BDC. 
Introducing the temperature ratio t = Tf,/T c and the minimal tem¬ 
perature ratio to > 1 , above which the engine is providing a power 
output, the gross power can be approximated with 

P *-“(tP/.T)(^r), (22) 

where a, in general, is a function of geometry parameters, working 
gas properties, the operating frequency /, filling pressure p/ and 
temperature T. 

The total internal power loss — also a function of /, p/, T, geom¬ 
etry parameters and working gas properties — can be written as 
P/,int = b(J,pf,T). The net power output of the engine is given through 

Pn=Pg~ Pi,in f (23) 

The influence of the temperature on a, and b is due to the 
temperature dependency of the working gas properties. Assuming 
a perfect gas, the net power becomes 

P n = a (/, p f ) - b (/, Pf ). (24) 

For a fixed frequency and filling pressure, the coefficients a, b 
and the minimal temperature ratio to can be obtained from the 
experiment. To do so, the net power output was measured for 
filling pressures of 6 bar, 8 bar and 10 bar, and frequencies from 
4 Hz to 12 Hz. For each pair of (f, p/) the engine was motored at 
heat input temperatures between 40 °C and 600 °C. The measured 
net power P n is plotted versus the temperature ratio t. For/= 8 Hz 
the results are shown in Fig. 12. For all investigated filling pres¬ 
sures the net power becomes zero at t = 1.62 ± 0.1. Since the 
internal power loss is a function of the filling pressure, and for all 
investigated pressures the net power equals zero at t = 1.62, the 
internal power loss has to be zero at i = 1.62, and thus 
P„(t = 1.62) = Pg(T = to) = 0, resulting in to = 1.62 =t 0.1. Hence, the 
engine needs a minimal temperature difference of T/ r -T c = 185 K 
to provide a power output. The net power as a function of t is 
adapted with Eq. (24), where for the minimal temperature ratio 
the value tq = 1.62 is used. 



Fig. 11. Calculated pV-diagram at a filling pressure of 10 bar and a operating frequency 
of 8 Hz. 



Temperature ratio 


Fig. 12. Fit of the measured net power output as a function of the temperature ratio for 
different filling pressures and a operating frequency of 8 Hz. 

With the fitted coefficients a and b, the total internal power loss 
Pi ,int and the gross power P g can be calculated. 

For filling pressures of 6 bar, 8 bar and 10 bar, and operating 
frequencies higher than 8 Hz, the internal power loss depending 
on the operating frequency is shown in Fig. 13. At each frequency, 
the values of the internal power loss are averaged and adapted 
with a function P/j nt f 5 . For operating frequencies lower than 
8 Hz the internal power loss is almost zero. Exceeding this fre¬ 
quency, it increases non-iinearly with the frequency, but is still 
quite low. 

For Th = 600 °C, and filling pressures between 6 bar and 10 bar, 
the gross power dependency on the operating frequency is shown 
in Fig. 14. At all pressures it is increasing linearly with the operating 
frequency. The gross power also depends linearly on the filling 
pressure. Thus, the coefficient a is directly proportional to the 
product of / and p/. This result is obvious since the gross power is 
proportional to the oscillating mass in the engine. The mass flux 
itself is directly proportional to the frequency and the filling 
pressure. 



Frequency [Hz] 

Fig. 13. Measured internal power loss as a function of the operating frequency for 
different filling pressures. To avoid too many lines in the plot the error bars are only 
shown for the averaged internal power loss. 
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Frequency [Hz] 

Fig. 14. Measured gross power output as a function of the operating frequency at a 
heat input temperature of 600 °C and different filling pressures. 

Fig. 15 shows the gross power output for a filling pressure of 
10 bar and heat input temperatures between 40 °C and 600 °C. In all 
cases, the observed gross power output increases (7/, > 200 °C), 
respectively decreases (Th < 200 °C), linearly with the operating 
frequency. The engine is not operating for heat input temperatures 
lower than 200 °C. 

8. Influence of design parameters on the engine performance 

In Ref. [8] it was shown that the engine performance strongly 
depends on the compression ratio and distribution ratio (see Sec¬ 
tion 4.2 for their definitions). In this work these theoretical pre¬ 
dictions are investigated experimentally and compared to the 
results of the simulation model. Additionally the influence of the 
cold heat exchanger area is studied experimentally. The model 
calibration (see Section 5) is kept at the values of the previously 
investigated engine geometry. Thus, the calculation results are 
quantitatively not comparable to the experimental observations. 
Nevertheless, the tendency of the performance by changing the 
engine geometry can be reproduced by the model as well. 



Frequency [Hz] 


Fig. 15. Measured gross power output as a function of the operating frequency at a 
filling pressure of 10 bar and different heat input temperatures. 



Frequency [Hz] 


Fig. 16. Measured net power output as a function of the operating frequency for 
different engine strokes. 

8.3. Compression ratio 

The engine was motored for strokes between 20 mm and 
60 mm. Fig. 16 shows the resulting measured net power outputs for 
operating frequencies between 4 Hz and 12 Hz. The highest values 
are observed for a stroke of 60 mm. With the experimental set-up, 
power outputs at higher strokes than this value were impossible to 
measure reliably. 

A comparison of the measured and calculated net power out¬ 
puts depending on the compression ratio according to Eq. (17) is 
given in Fig. 17. The operating frequency of the engine is 8 Hz, which 
is slightly lower than the net power’s peak frequency. Observation 
and calculation show an increasing net power for higher 
compression ratios. Thus, on a qualitative basis, the model re¬ 
produces well the observed engine behavior. For 7), = 400 °C the 
mismatch between calculation and observation is higher due to 
lower power outputs and thus a stronger influence of the modeling 
error. The experiment shows an optimal compression ratio of 0.47 
(50 mm stroke) for T/, = 400 °C and 0.56 (60 mm stroke) for 



Compression ratio 

Fig. 17. Calculated and measured net power output as functions of the compression 
ratio for different heat input temperatures. 
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Frequency [Hz] 

(a) T h = 200 °C 



Frequency [Hz] 

(b) T h = 400 °C 



Frequency [Hz] 

(c) T h = 600 °C 


Fig. 18. Measured net power output as a function of the operating frequency for 
different distribution ratios and heat input temperatures. 



Fig. 19. Calculated and measured net power output as functions of the distribution 
ratio for different heat input temperatures. 


Th = 600 °C. This observation is in accordance with the theoretical 
predictions in Ref. [8], Fig. 9. The compression ratio introduced in 
Ref. [8] is the compression ratio defined in Eq. (17) plus one. 
However, higher heat input temperatures allow higher compres¬ 
sion ratios and a higher power output. 

8.2. Distribution ratio 

To study the influence of the distribution ratio three re¬ 
generators with lengths of 111 mm, 122 mm and 143 mm are 
manufactured and placed into the 230 mm long stainless steel 
tube. The resulting lengths of the pulse tube are 76 mm, 68 mm 
and 44 mm. The flanges of the 122 mm regenerator version are in 
total 3 mm longer than in the two other regenerator versions. 
With these dimensions of pulse tube and regenerator the dis¬ 
tribution ratios according to Eq. (18) are rj = 1.2 (short regen¬ 
erator), r <i = 1.5 (medium regenerator) and r<j = 3.2 (long 
regenerator). For the motored engine, the measured net power 
outputs at heat input temperatures of 200 °C, 400 °C and 600 °C 



Frequency [Hz] 

Fig. 20. Measured net power output and efficiency as functions of the operating fre¬ 
quency for different cold heat exchanger areas. To avoid too many lines in the plot the 
error bars are only shown for the measured values at A chx = 112 cm 2 . 
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are shown in Fig. 18. For all temperatures the short regenerator 
always leads to the lowest power outputs. At Tj, = 400 °C the 
medium and the long regenerator lead to comparable power 
outputs. For lower temperatures the long regenerator becomes 
advantageous and for higher temperatures the medium regen¬ 
erator. An ideal engine having a high value of r<j (long regener¬ 
ator) has a high thermodynamic asymmetry and is according to 
[8], Fig. 10 functioning best. In the real engine the long regen¬ 
erator is accompanied by a short pulse tube. Thus, the hot heat 
exchanger is closer to the cold heat exchanger. A higher portion 
of the applied thermal power flows directly to cold heat 
exchanger and is released unused. This power loss increases for 
higher heat input temperatures and leads to the observed better 
engine performance with a medium regenerator length. 

In Fig. 19 the calculated and measured net power outputs 
depending on the distribution ratio for an operating frequency of 
8 Hz are compared. For all temperatures the calculations predict an 
even stronger dependency of the net power on the distribution 
ratio than the experiment shows. The best match is given for the 
medium size regenerator. The model was calibrated for this case. 
Nevertheless, the existence of an optimal distribution ratio is 
evident from measurement and calculation. 

8.3. Area of the cold heat exchanger 

It was claimed that the cold heat exchanger chxi needs a certain 
heat transfer rate for optimal engine operation. To check this its 
heat transfer area A C hx was modified experimentally. Three cold 
heat exchangers were manufactured with heat transfer areas of 
56 cm 2 , 112 cm 2 and 163 cm 2 , respectively. For T/, = 400 °C the 
measured net power output as well as the net efficiency are shown 
in Fig. 20. 

At frequencies lower than 8 Hz the influence of the cold heat 
exchanger area on the engine performance is small. Having higher 
operating frequencies the heat exchanger with an area of 112 cm 2 
delivers best performance results. The observation shows the 
importance of a careful design of the cold heat exchanger. In 
contrast to Stirling engines a maximized heat transfer capability 
does not lead to the best pulse tube engine performance. 

9. Conclusions 

The present study clarifies the working principle of the pulse 
tube engine and shows its potential and limitations. The developed 
numerical and experimental tools could be used successfully to 
built up a pressurized working engine and characterize its mode of 
operation under different heat input temperatures, filling pressures 
and operating frequencies. Model and experiment deliver compa¬ 
rable engine performances. 

The shaft power, the net power and the gross power of the 
engine could be obtained from measurement and calculation. 


Power losses were quantified. The investigated engine provides a 
maximum net power output of about 6 W at a frequency of 9 Hz and 
a net efficiency of about 8%. 

Theoretical predictions of a prior study concerning the influence 
of design parameters could be proven experimentally. The found 
optimal design parameters will help to develop future pulse tube 
engines. 

This work forms a basis for further studies of that kind of heat 
engine. The measured and calculated engine behavior for different 
heat input temperatures shows that the pulse tube engine has a 
potential application at low temperature heat sources (about 
300 °C). At this temperature level its practical performance will be 
comparable to other prime movers, such as Stirling engines. The 
theoretical disadvantage of the pulse tube engine process based on 
an intrinsically irreversible thermodynamic cycle could be 
compensated by its practical advantages, such as a comparably low 
flow friction power loss and the simple mechanical gear system 
associated with its low production costs. 
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